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We consider the problem of spin relaxation in a two-dimensional electron gas (2DEG) in a per- 
pendicular magnetic field. We assume that the spin relaxation is induced by the Rashba spin-orbit 
(SO) interaction, which appears due to the inversion asymmetry of the confining potential. Our 
solution is based on a microscopic evaluation of the spin density response function of the 2DEG with 
impurities and SO interaction. We derive explicit expressions for the transverse and longitudinal 
spin relaxation rates. Our analysis shows, in particular, that the spin relaxation rates exhibit mag- 
netoquantum oscillations, which are analogous to the Shubnikov-de Haas oscillations of the electrical 
resistivity. These oscillations can be observed, for example, in time-resolved optical experiments. 
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I. INTRODUCTION 

The study of spin relaxation in semiconductor 2DEGs 
is an important area of the emerging field of spintronicsX 
The dominant spin relaxation mechanism in these sys- 
tems can typically be associatedSi with the Rashba SO 
interaction;^ which exists due to the inversion asymmetry 
of the confining potential in semiconductor heterostruc- 
ture based 2DEG systems. A simple semiclassical pic- 
ture of the spin relaxation, due to Dyakonov and Perel 
(DP)ji2i4 is that since the Rashba interaction has the form 
of a momentum-dependent magnetic field Az • [r x p] , it 
induces precession of the spin of a moving electron. This 
precession leads to the randomization of the spin at long 
times. The DP spin relaxation rate is given, to leading 
order in A, by the simple expression I/t^^ ^ {Xpp/Tif'T, 
where pp is the Fermi momentum and r is the elastic 
scattering time. Quasiclassically, Xpp/fi is simply the 
spin-precession frequency, associated with the Rashba 
field. The distinguishing feature of the DP spin relax- 
ation is that the relaxation rate is smaller in more dis- 
ordered systems, since impurity scattering disrupts the 
Rashba spin precession by randomizing the electron's mo- 
mentum. 

It is well-known that transport processes in 2DEG 
systems are strongly influenced by the application of a 
perpendicular magnetic field. The perpendicular field 
quantizes electron's energy spectrum, which manifests 
spectacularly in Shubnikov-de Haas (SdH) oscillations 
of the resistivity and eventually leads to the quantum 
Hall effect. In this paper we address the question of 
how a perpendicular magnetic field influences the Rashba 
SO interaction-induced spin relaxation in semiconductor 
2DEG systems. This problem has been addressed before 
by several authors, both in the semiclassical^iSi^ and 
in the quantum limitSji^iiSiiiiiSiiii^ but a complete and 
rigorous analysis is still lacking. 

In this paper we will provide such an analysis. We will 
mainly be interested in the regime of moderate magnetic 
fields, such that uJcT ^ 1 (here ojc is the cyclotron fre- 
quency), and Ei? 3> hujc- Our analysis can, however, be 
generalized to the quantimr Hall limit, when ep ^ fiujc 



and oJcT is large. Our theory is based on a microscopic 
evaluation of the spin density response function of a dis- 
ordered 2DEG system with the Rashba SO interaction 
in a perpendicular magnetic field. We obtain explicit an- 
alytical expressions for the longitudinal and transverse 
spin relaxation rates: 
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which reduce to the well-known Dyakonov-KachorovskiiS 
expressions in the zero-field limit. We show, in particu- 
lar, that the application of a perpendicular magnetic field 
leads to quantum oscillations of the spin relaxation rate, 
which are analogous to the well-known quantum oscil- 
lations of transport coefficients. These oscillations arise 
from the magnetic field dependence of the scattering time 
T, which will be calculated below. 

The paper is organized as follows. In section^lwe cal- 
culate the disorder-averaged Green's function of a 2DEG 
system with the Rashba SO interaction in a perpendicu- 
lar magnetic field using the self-consistent Born approxi- 
mation (SCBA). In section iHll the spin-density response 
function of our system is calculated by summing SCBA 
self-energy and ladder vertex corrections to the bare po- 
larization diagram. The transverse and longitudinal spin 
relaxation rates are then extracted from the poles of this 
response function. 



II. DENSITY OF STATES AND DISORDER 
SELF-ENERGY 

The simplest single-particle Hamiltonian describing 
the dynamics of electrons in a semiconductor 2DEG in a 
perpendicular magnetic field can be written as: 



= h Az • 

2 m 



T X TT 



(1) 



where tt — p + (e/c)A is the kinetic momentum, r is the 
spin operator and = gfJ'sB is the Zeeman energy. The 
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second term in Eq.(^ describes the Rashba SO interac- 
tion. Hamiltonian Hq can be easily diagonalizedi^iiSiii 
if one notices that the Rashba term mixes only neigh- 
boring Landau levels with opposite spin directions. The 
eigenstates are thus given by: 



Green's function. This is a very good approximation for 
the high Landau level filling factors that we are assum- 
ing. SCBA equation for the retarded disorder self-energy 
in our case reads: 



\n,a) = Una\n,i) + Vna\n - 1,T), 



(2) 



^na(^) ~ ^ ] ^na.n'a'Gn'a'i^)j 



(8) 



where a = 1, 2 and |n, a) is the n-th Landau level eigen- 
state with spin a =t, i- The corresponding eigen-energies 
are: 



where 



k 



e„Q = huJcn + (-1)" 




and 



2X'^mhuJcn, (3) 



(10) 



where ujc = eB/mc. Note that n = 1, 2, ... for a — I, but 
n = 0, 1, ... for a = 2. The amplitudes Una and w„a are 
given by: 




is the retarded SCBA Green's function. Angular brackets 
in Eq.lQ denote disorder average. 

Matrix elements of F can be easily evaluated for 6- 
function impurity potential and are given by: 



A^y/i + 2\^mfLUJcn 



27rTo 



(11) 



(-l)°+i(^c^c + A,)/2 
2 ' y/{hujc + A^)2/4 + 2\^mhuOcn ' 

(4) 

Using the above basis of single-particle eigenstates, the 
Hamiltonian of the 2DEG system in the presence of im- 
purity scattering potential can be written as: 



where tq is the elastic impurity scattering time in the 
absence of the magnetic field. 

Eq.ljSJ is too complicated to be solved analytically 
without approximations. Fortunately, only very minor 
simplification is needed to make this equation solvable. 
The main complication is the dependence of the matrix 
F on the Landau level index n. This dependence is, how- 
ever, inessential, since it appears entirely due to the term 



u — c „ I T/ , ,(h h'\n^ n 2A2m?iWc'T. in Eq.EJ. The dependence on n in this term 
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nak.n' a' k' 



(5) 

where index k denotes orbit-center quantum numbers in 
the Landau gauge and 

Ka,„'a'(fc,fc')=E / rfrv|/;,,(m)T/(r)vE'„,,,fc,(m), (6) 

are the matrix elements of the impurity potential V{y) 
in the basis The wavefunctions 4'nafc(ro') are given 
by: 



can be safely ignored, since the Rashba SO interaction 
can be assumed to be a weak perturbation, in the sense 
that Xpp ep. In this case we can make a replacement 
n —>■ ep /hiUc in the Rashba term in Eqs. Q and Q , which 
makes the amplitudes Una and Vna and, consequently, F, 
independent of n. This simple approximation allows us 
to solve Eq.® analytically. 

It is convenient to introduce the following new nota- 
tion. Let 



A 



*na/c(r T) = W„a0n-l,fc(r), 
*na/c(r i) = Una(l>n,k{r), 




2X^mep. 



(12) 



(7) 



where 4>n,k{^) are the Landau- level eigenfunctions in the 
Landau gauge. 

Our analysis of the spin relaxation in a system, de- 
scribed by Eq.Q, is based on a calculation of the spin 
density response function, which is similar to Ando's cal- 
culation of the conductivity of a 2DEG in a perpendicular 
field. Our calculation is technically more complicated 
than Ando's due to the presence of the SO interactions 
and also due to the fact that vertex corrections to the 
polarization bubble, which vanish in the case of the con- 
ductivity calculation, are crucial in our case. 

Following Ando, we use the self-consistent Born 
approximation (SCBA) to find the disorder-averaged 



Then e„a = T^^cn + (-l)''A/2. Also define: 
ui — icos(i9), vi — sin(i9), 
U2 — isin(i9), V2 = — cos(i?), (13) 

where 



cos(?9) 



1 hujc + Az 



2 2A 



(14) 



Using this notation, the matrix elements of F are given 
by: 



„ Tl UJr f It,. 

Fn =F.. = ^h--sin^(2.) 



Fl2 — F21 



h^UJr 



TTTo 



sm'^{2d). 



(15) 
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Let us also make the usual assumption that the real part 
of the disorder self-energy can be absorbed into the chem- 
ical potential. Then SCBA equation can be written as: 



where 



(?a(e) 



:iniG,l(e), 



(16) 
(17) 



is the density of states and I = ^Jhc/eB is the magnetic 
length. Writing out Ea. (|17|l explicitly, we have: 
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E 



1 



ImS^(e) 



[e - huo.n - (-l)-A - hLU,5a,i]- + [ImE^(e)]" 

(18) 

The sum over Landau level indices in the above equation 
can be done using Poisson summation formula*^^ 



n=0 



/'OO 



(19) 



Up to this point, our calculations were applicable to mag- 
netic fields of general strength. To proceed further, we 
will restrict ourselves to the regime when epT/h ^ 1 
and ep ^ hujc- (The calculation we present can be done 
in other regimes as well, but the actual procedure will 
be a little different.) This allows us to neglect the term 
/(0)/2 in Ea. (|19|l . Then, changing the integration vari- 
able to a; = hujcTi — t+ (— 1)°A — fnjJc5a,i and extending 
the lower limit of integration to —oo, which is justified 
by the above assumptions, we obtain: 



Qa{e) 



1 



27r2 



dx 

TlbJr 



X exp ( 2'!Tip 



P— — 00 

x + e- (-1)"A - hujcSa.i 

TlLJr 



Calculating the above integral, we obtain 
ea(e) 



P«(e)]'' 
(20) 



2-Kn 

X cos I 2ttp 



1 + 2 ^ g-^P/'^cTo 



htOr 



(21) 

where we have replaced ImE^(e), which appears in the 
exponent after integration, by its zero- field value —h/2TQ. 
The oscillatory term in the above expression is slightly 



different from the usual expression, which is periodic in 
hujc, due to the presence of the Rashba and Zeeman split- 
tings. In what follows we will ignore this modification of 
the density of states, putting A = and = in the 
oscillatory term. This does not mean, of course, that we 
are ignoring SO interactions altogether, since they still 
enter in the matrix elements of F. Then the density of 
states becomes independent of the index a and is given 
by the following simple expression: 



go 



np 



(22) 



where gii(e) = P2(e) = f?(e)/2 and go = m/h^ir is the 
zero-field density of states. If the magnetic field is not 
too large, i.e. if uJcT < 1, only the first term in the sum 
over p in Ea. 12211 needs to be retained. Then the oscil- 
latory term in the density of states is purely sinusoidal. 
Substituting Eq.(|21|) in F,q.^ we find that SCBA self- 
energy is also independent of the index a and is given by: 



ImS^(e) 



2r(6) 



2mTQ 



(23) 



The above disorder self-energy reduces to the usual ex- 
pression —h/2To in the zero-field limit. At finite fields, 
the oscillatory term in the density of states in Ea. H23() 
is at the origin of a number of magnetoquantum oscil- 
lation phenomena, including the SdH oscillations of the 
electrical resistivity. As will be shown later, it also leads 
to similar magnetoquantum oscillation effects in the spin 
relaxation rate. 



III. SPIN DENSITY RESPONSE FUNCTION 

We can now calculate spin density response function of 
a disordered 2DEG system with the Rashba SO interac- 
tions in a perpendicular magnetic field. Transverse and 
longitudinal spin relaxation rates can be found from the 
poles of the corresponding spin density response func- 
tions. The calculation is similar to the analogous cal- 
culation that can be done at zero fieldiSS We will use 
h = 1 units in the rest of this section, except for the final 
results. 

We start from the general density matrix response 
function: 

Xo'i<T2,CT3(T4 (r — r,t t ) — 
- teit-t'){[gl,^ir,t),g,.^^,^ir',t')]), (24) 

where pcriCT2(r,t) — '^l,^{r,t)^cri{^,t) is the generalized 
density operator, whose expectation value is the density 
matrix. As usual, the calculation is most conveniently 
done by analytical continuation of the imaginary time 
response function. Fourier transformed imaginary time 
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(a) 






(b) 



FIG. 1: (a) SCBA equation for the disorder-averaged Green's 
function; (b) Equation for the vertex function D. Dashed 
lines denote impurity potential correlator {V{r)V{r')). 



response function can be written as: 

XcriO-2,(T3 £74 

where 



(25) 



0'1<T2 .<T3 0-4 



(r — r', ioj, iuj + iH.) 



{r,r' ,iuj + ift)) (26) 

is the polarization bubble diagram and G are imaginary 
time Green's functions. P can be calculated by sum- 
ming all the SCBA diagrams for the disorder self-energy 
and all the ladder vertex corrections to the polarization 
bubble, which constitutes a conserving approximation for 
the density matrix response function. The result can be 
written in matrix notation as P = P^D, where is 
the "bare" polarization bubble with only the self-energy 
corrections included: 



Qa3<Ji (r', r, iL^)Qa2ai (r, + i^), 



(27) 



with Q here being the disorder-averaged Green's func- 
tion. The vertex part, or, as it is often called, diffusion 
propagator, D satisfies the following matrix equation: 

P> = 1+7P"D, (28) 

where 7 = (F^(r)}. The solution of this equation is 



D= [I-7P"]"'. 



(29) 



We can now calculate the retarded real time response 
function by analytically continuing Ea. l|25|) to real fre- 
quencies. The result can be written as an integral along 
the branch cut of the SCBA Green's function: 

t°° de 

X(q,ri)- / —nF{€)[P{q,e + iTj,e + n + t7j) 
J-oc 27ri 

— P(q, e ~ iri,e + n + irj) + P(q, e fl — iri,e + iif) 

- P(q,e-r!-i77,e-z77)]. (30) 



At low frequencies and low temperatures Ea. H30() can be 
simplified to: 

X(q, n) = gp(q, ep -i7j,eF + n + i^) + (31) 

Thus, the problem reduces to calculating the matrix 
P''(q, ep — ir^^ep +^ + i'n)- We will in fact be interested 
only in the relaxation of uniform spin polarization and 
thus will put q = henceforth. The calculation of P° 
in our case closely resembles the analogous calculation in 
the zero-field problem. 

We start from the definition of P° Ea. (|27|l . using the 
disorder-averaged Green's functions calculated in section 

m 



1 

L 



<Ti<T2,o-3<T4 (^F - i?7, CF + + +ir]) 

ri n h Ti ' n' 1< 



nak,n' a' k 



i/2T 



(32) 



where r = T{ep) and the wavefunctions are given by 
Eq.(I7|). The spatial integrals in the above equation can be 
easily evaluated using the orthonormality of the Landau 
gauge eigenfuctions and are given by: 



dr4':„fe(r,i)*„,„,fc,(r,i) 
«fr*Lfe(r,T)vI/„,,,,,(r,i) 
dr*:,fe(r,i)«'„,,,fc,(r,T) 



U*aU^,Sn,n'Sk,k' , 
ulv^,6n,n' -lSk,k' -{33) 



The sums over Landau level indices in Ea. (|32|l can be 
done as follows. We first decouple the product of two 
Green's functions appearing in H32I) as: 



1 



{ep - ena - i/2T){ep + VL- e^'a' + ij^r) 
1 

VL + (-1)"A - (-1)"' A -t- ijr 

1 1 



(34) 



The real parts of the decoupled sums over Landau level 
indices can then be neglected by our assumption that 
tp ^ hiOci ^PFi The imaginary parts are simply pro- 
portional to the density of states at Fermi energy Q{ep). 
Thus we obtain the following expressions for the matrix 
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elements of P": 



rates are given by: 



sin2(2?9) 



[2/o(r!)-/+(r!)-/_(r!)] 



7^ 



/o(r! + c^e)+sin^W/+(r! + c^c) 

+ cos^(i?)/_(r! + c^,), 
sin2(2?9) 



Ti.Ti - 2 

4 



/o(r!-c^e)+sin4(i?)/_(r!-w,) 
+ cos\d)f+in-uJc), (35) 



iT,iT ~ 2 

A 



where functions /o and j± are given by: 

1 



1 



1 - iflT ± I At' 



(36) 



All other matrix elements of P'^ vanish. 

To calculate the transverse and longitudinal spin relax- 
ation rates we need to find the poles of the spin density 
response function, or, equivalently, zeroth of the inverse 
diffusion propagator D~^. It is most easily done by trans- 
forming D^^ to the physical space of charge and spin 
density components, which is accomplished by multiply- 
ing it on both sides with Pauli matrices: 



(37) 



where [3 — c, x, y, z and r"^ is the identity matrix. Thus 
we obtain: 



D-i(f^) = l-/o(f^) 



sin2(2?9) 



D+'Lin) = 1- 



sin2(2i?) 



fo{n + u;,)-sin^{i9)f+{n + Lo,) 



cos''(t?)/_(f7 + Wc). 



(38) 



Clearly D^^{0) — 0, which is a consequence of charge 
conservation. For the spin part of the inverse diffusion 
propagator, to leading order in Xppr/h and fir, and as- 
suming that Az <S hujc, we obtain: 

Tz 

T^^Dilifi) = -in + in^^Az + —, (39) 

where the longitudinal and transverse spin relaxation 



1 

Tz 
1 



1 + (c^cT)2 
1 

277' 



(40) 



Ea. (|40|) is our main result. The form of these ex- 
pressions is rather similar to the well-known Dyakonov- 
Kachorovskii expressions^ for the spin relaxation rates in 
the zero field case, except for the factor 1/[1 -|- {lUct)^], 
which describes the suppression of spin relaxation by the 
perpendicular field. The mechanism of this suppression 
can be understood in quasiclassical terms by imagining 
electrons moving along classical cyclotron orbits and un- 
dergoing impurity scattering and the Rashba spin pre- 
cession. Orbital motion affects the spin precession, since 
the direction of the Rashba field is always transverse 
to the direction of electron's velocity, and thus changes 
as the electron moves along a circular cyclotron orbit. 
However, if cOcT ^ 1, impurity scattering randomizes 
the direction of electron's motion and the influence of 
the magnetic field on the orbital motion is then negligi- 
ble. On the other hand, when uJcT 3> 1, electrons can 
move around cyclotron orbits almost freely and thus the 
Rashba spin precession is averaged out. This picture of 
the suppression of the DP spin relaxation by a perpen- 
dicular magnetic field has been discussed before by sev- 
eral authorsi^iSii However, we are not aware of a rigorous 
derivation of a simple analytical expression for the spin 
relaxation rate, such as Ea. H4U|) . 

Probably the most noteworthy feature of our result is 
the oscillatory dependence on the magnetic field, which 
appears due to the oscillatory term in the scattering rate 
1/t in Eas. (|22|) and (|23|l . If uJcT < 1, we may neglect 
the (uJct)'^ term in the denominator in Ea. H4U|l and leave 
only the first term in the sum over p in Eq. (|22|l . In this 
case spin relaxation rates will exhibit purely sinusoidal 
SdH-like oscillations. Similar effects have been discussed 
earlier in the context of nuclear spin-lattice relaxation 
in quantum Hall systems However, the possibility 
of such magnctoquantum oscillations in the SO-induced 
spin relaxation rates has never been discussed before. 
Such oscillations have in fact been observed in recent ex- 
periments on spin relaxation in InGaAs quantum wellsiS^ 

In conclusion, we have considered the problem of spin 
relaxation in a 2DEG with the Rashba SO interactions in 
a perpendicular magnetic field. We have given a rigorous 
microscopic derivation of the transverse and longitudinal 
spin relaxation rates from the poles of the spin density 
response function in the limit S> tkuc- It is possi- 
ble to extend our calculation to other regimes, including 
the quantum Hall regime ep ^ fiujc- Our results demon- 
strate that magnetoquantum oscillation effects should be 
observed in spin relaxation phenomena in 2DEGs. One 
of the advantages of our approach is that it can be easily 
extended to study another very interesting and still un- 
explored question of what is the effect of Coulomb inter- 



6 



actions on the spin relaxation in semiconductor 2DEGs. 
This question wiU be addressed in a future pubUcation. 
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